Ab initio study of misfit dislocations at the SiC/Si(001) interface 
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The high lattice mismatched SiC/Si(001) interface was investigated by means of combined classical 
and ab initio molecular dynamics. Among the several configurations analyzed, a dislocation network 
pinned at the interface was found to be the most efficient mechanism for strain relief. A detailed 
description of the dislocation core, and the related electronic properties, are discussed for the most 
stable geometry: we found interface states localized in the gap that may be source of failure of 
electronic devices. 
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In recent years the development of several sophisti- 
cated epitaxial growth techniques has allowed to produce 
new interesting materials and heterostructures, even in 
the presence of high lattice mismatch between the con- 
stituent phases. Furthermore, high-resolution electron 
microscopy (HREM) and other similar techniques deep- 
ened the characterization of such systems, providing a 
description of the interface quality and epitaxial relations 
of the sublattices. Nevertheless, the atomic structure and 
the chemical environment at the interface, which deeply 
affect its physical properties, still remain hardly accessi- 
ble to any experimental technique. 

Atomistic simulations represent a powerful tool to 
complement the experimental data. Though recent theo- 
retical works appeared on the study of lattice mismatched 
interfaces for some specific system involving ceramic ma- 
terials the ab initio study of a mismatched semicon- 
ductor/semiconductor interface, where the complex cova- 
lent bonding requires an accurate description of charge 
transfer is far to be complete @. Among all the 
possible high misfit heterostructures, we focussed on the 
(001) interface between cubic silicon carbide (/3-SiC) and 
silicon (Si) without lacking of generality. This system, 
characterized by a peculiar ^20% difference in lattice 
parameters between the two constituents, can be con- 
sidered as a template of high lattice mismatched hetero- 
junctions, where an ab initio approach is still affordable. 
Furthermore, it is technologically interesting, because of 
the outstanding physical properties of /3-SiC, a possible 
candidate for applications in harsh environment devices 
. Some recent results on the study of this heterostruc- 
ture ^ evidenced that the interface can be prepared so 
to be locally abrupt, with a network of pure edge mis- 
fit dislocations localized at the interface to accommodate 
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the extreme lattice mismatch. The formation of a dislo- 
cation network as an efficient mechanism of strain release 
usually occurs in heterostructures with a mismatch larger 
than -10% in. 

The occurrence of edge dislocations can be treated in 
the frame of the well established near coincidence lat- 
tice model 8], which asserts that the lowest interfacial 
energy configurations are obtained when a perfect coin- 
cidence site between two dissimilar structures is realized. 
Given two materials with lattice parameters ai and a2, 
a perfect coincidence site occurs when ai/a2=m/n, with 
m and n positive integers. For the /3-SiC/Si(001) sys- 
tem, m=5 and n=4 respectively: since m=n-|-l, only one 
extra lattice plane is needed along each of the two prim- 
itive directions of the coincidence lattice cell to release 
almost all the misfit strain, in agreement with HREM 
images [1, Q . A network of misfit dislocations will thus 
be generated in the SiC film starting at the interface. 

In order to characterize the heterostructure and pro- 
vide an accurate description of its structural and elec- 
tronic properties, we performed a combined classical and 
ab initio molecular dynamics study of different disloca- 
tion core structures: the minimum energy system, i.e. 
the most efficient way for strain release, is obtained for 
dislocations pinned at the SiC/Si interface, originating 
at a sub-stoichiometric, C terminated layer. 

To reproduce the experimental situation, we consid- 
ered a periodically repeated multilayer: the interface 
was built by matching a p(4x4)-Si(001) with a p(5x5)- 
SiC(OOl) slab, at the Si lattice parameter optimized for 
bulk calculations with the chosen potentials (either clas- 
sical or ab-initio) . This model contains two edge disloca- 
tions in the supercell, one along each of the two directions 
perpendicular to the interface, and the periodicity exper- 
imentally observed is thus respected. In Fig. [T] the inter- 
face structure is schematically represented and the coin- 
cidence lattice supercell, which we considered in our cal- 
culations, is evidenced. At the interface, the dislocation 
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directions are not equivalent, because of the zincblende 
stacking: this renders inequivalent the dislocation cores. 

We have performed a structural optimization for sev- 
eral possible interface configurations, varying both geom- 
etry and stoichiometry, via classical and ab initio Molec- 
ular Dynamics (MD) methods in order to get the most 
stable structure. The relative energy of systems contain- 
ing different numbers of atoms has been discussed in the 
grand canonical scheme Q. 

Classical MD for several interface configurations al- 
lowed us to select the core dislocation structures and slab 
size for the ab initio simulations. The Si-Si, Si-C and C- 
C interactions were modelled by means of the empirical 
Tersoff's potentials [10]. It has recently been demon- 
strated that these potentials give formation energies and 
properties of native defects in /3-SiC(001) in good agree- 
ment with ab initio calculations [ll[ . We performed tests 
varying the stoichiometry at the interface and the core 
dislocation distance from the interface. In order to pre- 
vent the interaction between the interface and the slab 
surfaces, we considered thick slabs (36 SiC layers on top 
of 36 Si layers: in the following, we will indicate the sys- 
tems as N/M, where N (M) stems for the number of SiC 
(Si) layers, thus the slab described above will be labelled 
as 36/36). All the atoms were allowed to move, except 
those belonging to the outermost layers, which were con- 
strained to relax in a collective way. Our results indicate 
that configurations with C atoms at the interface are al- 
ways more stable. A pseudomorphic SiC layer on top 
of Si, before the dislocations start, is unfavorable. In 
agreement with experimental results , the favored core 
structures are those with the dislocations pinned at the 
interface p^]. Among all the considered geometries, we 
focussed on the most meaningful ones: the stoichiometric 
interface and the one which is lowest in energy. The latter 
is obtained from the stoichiometric solution by removal 
of all the C atoms along the two core dislocation lines (13 
atoms per supercell): this configuration gives an energy 
gain AE of 1.12 J/m^ (16.43 cV/ccll-area) and 0.48 J/m^ 
(7.07 eV/cell-area) for C-poor and C-rich conditions re- 
spectively in classical MD and for a 36/36 slab. 

In order to determine the minimum size required to 
get accurate results via ab initio techniques, we analyzed 
how the relative energies and atomic structures of these 
two systems vary as a function of the slab thickness, in 
classical MD calculations. We observed that the relative 
atomic distortion is less than 2 % after the 4th layer, 
when considering a 5/5 or a 36/36 layer slab; further- 
more, the value obtained for AE for the two configura- 
tions described above in a 7/7 slab reproduces correctly 
the asymptotic value of an infinite (36/36) system. When 
a 5/5 slab is considered, AE decreases by ~ 20% in C- 
poor conditions (i.e. AE = 12.75 eV/cell-area), due to 
a stronger surface-interface interaction obtained in the 
non-stoichiometric structure. These results, that provide 
a good estimate for the system stability, allowed us to 
consider only the 5/5 slab in ab initio calculations. 

The ab initio MD simulations were performed in 



the frame of density functional theory (DPT), in the lo- 
cal density approximation (LDA) [IJ]. The electronic 
wavefunctions (charge density) were expanded in plane 
waves with an energy cutoff of 40 (160) Ry; only the 
r point was included for the integration in the supercell 
Brillouin Zone (BZ). The electron- ion interaction was de- 
scribed by fully non local pseudopotentials [Tsl. Surface 
atoms at both sides were saturated with hydrogen, at a 
distance optimized via preliminary surface calculations 
All the atoms of the slab were allowed to move ex- 
cept the hydrogens and the Si (C) in the outermost layers, 
which were relaxed along the direction perpendicular to 
the surface only. In a subsequent test, this constrain has 
been released in order to evaluate the interaction between 
surfaces and interfaces. Structures were considered con- 
verged when forces acting on atoms were less than 10^** 
a.u. (0.005 eV/A) and energy varied by 10 ^ eV/atom. 

We observed that the classical simulations give good 
initial guess geometries for ab initio calculations, though 
lacking of accuracy in predicting structures where dan- 
gling bonds are present. In particular, at variance with 
the empirical simulation, the stoichiometric system re- 
constructs around the two dislocation cores, where dan- 
gling bonds are localized, if relaxed by the ab initio tech- 
nique. The core structures are thus characterized by 5- 
and 7-membered rings with C-C and Si-Si dimers at the 
SiC and Si interface layers respectively [lB|. The same 
energy ordering was however found in the two different 
methods: the sub-stoichiometric solution is the most sta- 
ble in the full physical range of variability of the C chemi- 
cal potential. The ab initio energy gain for this geometry 
is AE = 14.9 (5.6) eV/cell-area, in C-poor (C-rich) con- 
ditions. 

The relaxed structure for the lowest energy interface is 
reported in Fig. O The interplanar distance at the inter- 
face is lower (-13%) than the one obtained in bulk SiC: 
the deformation optimizes the SiC cell volume and bond 
distances at this high lattice mismatch. We observed 
that both Si-Si (in Si) and Si-C bonds (in SiC and at the 
SiC/Si interface) are stretched at the dislocation cores, 
compared to the bulk bond lengths, being elongated up to 
2.5 A and 2.0 A, respectively. The tetrahedral configura- 
tion is distorted too, with some Si-C-Si bond angles close 
to 120°. The bulk structure is quickly recovered when 
moving away from the interface and at the third atomic 
layer the residual distortions are already small (±3% in 
bond distances, interplanar spacing and angles). The de- 
formation induced by the interface and the dislocation 
cores has been also analyzed in terms of warping of the 
atomic layers (see right panel in Fig. [J). We found that 
the major deformations are localized in Si, since it has 
smaller elastic constants than SiC. The layer puckering 
decreases when moving aside the interface: this is in nice 
agreement with some recent experimental results |l7| on 
the structural characterization of SiC films grown on a 
Si(OOl) substrate, which evidenced an internal roughness 
of individual SiC planes fading away from the Si sub- 
strate. As noticed in Fig. [21 the sub-stoichiometric struc- 
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ture does not present any dangling bond: the Si atoms in 
the second SiC layer, which would have an unsaturated 
bond, reconstruct and form dimers 2.4-2.5 A long; this is 
also evidenced by the analysis of the electronic structure 
(see below) which proves the presence of bonding charge 
density localized between interfacial Si atoms. The re- 
moval of C atoms and stretched C dimers in the core 
allows to better accommodate the misfit strain, although 
maintaining an abrupt SiC/Si interface: this stabilizes 
the system with respect to the stoichiometric one. 

In order to have a qualitative indication of the inter- 
face/surface interaction, we performed further calcula- 
tions on the 5/5 slab, relaxing the surface atoms at both 
sides. The results of these simulations indicate that al- 
ready for the 5/5 system the interface is well described, 
and the interaction with the surface is small: the bond 
lengths around the dislocation cores change less than 3% 
and the energy difference between the stoichiometric and 
sub-stoichiometric slabs is modified by 18% for C-poor 
conditions. The major distortions obviously occurs at 
surface layers, that bent to accommodate the residual 
elastic strain (Az — 0.9 A). 

The interface structure determines the electronic prop- 
erties of the heterojunction. In particular, the presence 
of defects such as misfit dislocations can induce interface 
states, also called interface trapped charges, localized in 
the band gap, that may be responsible for device failure. 
We have analyzed the electronic structure for the most 
stable dislocation network at the /3-SiC/Si(001) interface. 
Although the number of dangling bonds is minimized in 
this system, the Highest Occupied (HO) and Lowest Un- 
occupied (LU) states are localized in the core of the edge 
dislocations, as a result of the large difference in charge 
transfer between Si-C and Si-Si bonds. These states lay 
in the interface band-gap, at 0.8 and 1.3 eV above the 
valence band top at F [l8|, thus they may be a pertur- 
bation to the electron mobility and optical properties of 
the SiC/Si heterostructure. In Fig. O the charge density 
plot of the HO state along the two dislocation directions 
is represented. The charge density is mainly localized 
on atoms of the [110] core dislocation, while no density 
is observed around the core in the perpendicular direc- 
tion. The opposite situation is found in other bonding 
states localized in the system forbidden gap and for the 
LU state (not shown) which charge densities pertain to 
the dislocation core existing along the [110] direction. 

In conclusions, using classical and first-principles cal- 
culations, we have presented a comparative analysis of 
the structural and electronic properties of edge disloca- 
tion networks at a high lattice mismatched semiconduc- 
tor/semiconductor interface. For the studied case, /3- 
SiC/Si(001), the core dislocations which were found as 
most stable, and which are thus responsible for strain 
relief in the heterostructure, are pinned at the interface, 
and occur at sub-stoichiometric C layers. The presence 
of these extended defects is responsible for occupied and 
unoccupied states, laying in the forbidden gap, which can 
severely affect the electronic properties of the system. 
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FIG. 1: Schematic view of the slab (left) and atomic positions 
of the unrelaxed C-Si interface projected along (001) (right). 
The computational unit cell (cubic coincidence lattice) is in- 
dicated in continuous lines and the edge dislocation lines [110] 
and [110] are evidenced. Black (grey) spheres represent C (Si) 
atoms. Only two planes at the interface are shown. 
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FIG. 2: (110) (left) and (110) (center) projection of the most 
stable interface structure. A projected view of the interface 
atoms is presented as an insert (right) to enhance the warping. 
Black (grey) spheres represent C (Si) species. Surface atoms 
are saturated with hydrogen (small white spheres). 



FIG. 3: Charge density contour plot of the highest occupied 
state projected along (110) (left) and (110) (right). 
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